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INTRODUCTION 


In  the  areas  of  ballistic  missile  defense,  communications,  electronic 
warfare,  and  missile  design  and  testing,  the  investigation  of  nuclear 
effects  on  electromagnetic  wave  propagation  is  necessary.  Predictive 
models  of  atmospheric  modification  in  the  10  to  100  km  altitude  region 
resulting  from  nuclear  events  require  knowledge  of  prenuclear  atmospheric 
composition.  It  is  within  the  10  to  100  km  altitude  interval  that  least 
information  is  available  concerning  atmospheric  composition. 

The  existence  of  a functioning  computer  program  to  model  atmospheric 
composition  [1]  makes  it  possible  to  provide  theoretical  support  to  ASL 
researchers  in  their  ongoing  efforts  to  study  the  conditions  and  composi- 
tion of  the  atmosphere.  A previous  report  [2]  was  released  to  present 
part  of  the  results  obtained  for  the  32°  N September  stratosphere.  Hence, 
as  a followup,  the  computerized  model  was  recently  utilized,  to  generate 
comparative  results  for  the  65°  N September  stratosphere. 

The  basic  chemical  reactions  used  are  the  same  as  partially  described  in 
[2]  and  more  completely  in  [3].  Photoreaction  coefficients  were  computed 
by  J.  L.  Collins  and  are  described  in  [4]. 

In  attempting  to  deal  with  transport  the  same  approach  as  discussed  in  [2] 
was  used  to  develop  a model  of  the  chemical  composition  in  the  September 
stratosphere  at  65°  N latitude. 

The  required  temperature  profile  was  obtained  from  R.  0.  Olsen  of  the 
Atmospheric  Sciences  Laboratory.  A literature  search  was  performed  by 
J.  L.  Collins,  formerly  of  the  University  of  Texas  at  El  Paso,  to  obtain 
initial  values  of  concentrations  for  some  key  species. 


COMPUTATIONAL  RESULTS 

A similar  format  for  presentation  as  was  used  to  describe  the  results  of 
the  32°  N September  stratosphere  is  followed  to  give  the  results  for  the 
65°  N September  stratosphere. 

Consequently,  the  categories  of  calculated  results  presented  consist  of: 
(1)  particle  densities  of  selected  important  species,  (2)  mixing  ratios 
by  volume  of  these  same  species,  and  (3)  reaction  rates  of  several  sets 
of  some  key  photolytic  and  chemical  reactions.  All  results  shown  are 
given  as  a function  of  altitude  between  10  and  50  km  at  noontime.  Com- 
putations were  done  at  5 km  intervals. 

From  the  nitrogen/oxygen  family.  Figs,  la  through  lk  are  computed  altitude 
profiles  of  the  particle  densities  of  ground-state  oxygen  atoms  (0(3P)), 
excited  oxygen  molecules  (02(1a)),  ozone  (03),  nitric  oxide  (NO),  nitrogen 
dioxide  (N02),  N03  radical,  dinitrogen  pentoxide  (N205),  nitric  acid 
(HN03),  nitrous  oxide  (N20),  along  with  nitrous  acid  (HN02),  and  electron- 
ically excited  oxygen  atoms  0(1D). 
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From  the  oxygen/hydrogen  family.  Figs.  U through  lo  are  the  computed 
particle  densities  for  water  (H20),  hydroxyl  radical  (OH),  hydrogen 
peroxide  (H2O2)  and  perhydroxyl  radical  (H02). 


Methane  (CH4),  formaldehyde  (HCHO),  »nd  carbon  monoxide  (CO)  from  the 
carbon/hydrogen/oxygen  family  of  reactions  are  shown  in  Figs.  Ip,  lq, 
and  lr. 

Existence  of  some  of  the  species  in  the  atmosphere  is  speculative;  but 
since  their  presence  is  required  by  photolytic  and  chemical  considerations, 
these  species  are  also  presented  here.  For  some  species  which  have  been 
detected  in  the  atmosphere,  the  measured  results  are  plotted  along  with 
the  computed  values  to  offer  a comparison  between  theoretical  calculations 
and  experimental  measurements. 

Mixing  ratios  for  the  same  data  as  la  through  lr  are  presented  in  Figs. 

2a  through  2r  to  allow  direct  consideration  of  the  fraction  of  the  atmo- 
sphere that  a given  constituent  represents.  This  ratio  of  the  numbers  of 
the  particles  is  termed  the  volume  mixing  ratio  and  is  computed  by  divid- 
ing the  individual  particle  density  by  the  total  particle  density  of  all 
the  species  present  in  the  volume  considered.  These  figures  express  the 
fact  that  species  which  play  dominant  roles  in  many  atmospheric  processes 
exist  in  only  extremely  small  proportions  in  comparison  to  other  less 
chemically  important  species. 

To  better  understand  the  role  of  various  chemical  reactions,  one  should 
consider  the  rates  at  which  given  reactions  occur  under  the  prescribed 
conditions.  Hence  Figs.  3a  through  3n  are  the  calculated  altitude  depen- 
dent rates  of  a number  of  reactions  for  the  65°  N September  situation. 
Although  these  rates  are  labeled  "reaction  rates,”  they  are  more  precisely 
defined  "volume  production  rates." 

A partial  list  of  the  reactions  is  shown  in  Table  2.  On  the  graphs,  some 
of  the  reactions  are  grouped  to  allow  immediate  observation  of  certain 
critical  relationships  between  species.  These  relationships  are  described 
in  [2]. 


THE  ATMOSPHERIC  CHEMICAL  KINETICS  MODEL 

The  design  and  calculational  method  used  in  the  computer  simulation  of 
stratospheric  composition  as  well  as  the  structure  of  the  34  photodisso- 
ciative  processes  and  115  chemical  reactions  are  discussed  in  [2]  and 
[3],  The  atmospheric  molecules,  atoms,  and  free  radicals  used  in  the 
modeling  effort  are  shown  in  Table  1.  A selected  list  of  reactions  with 
rate  constants  is  given  in  Table  2. 


TABLE  1 


ATMOSPHERIC  CHEMICAL  SPECIES  CONSIDERED 
IN  COMPUTATIONAL  MODEL 

0 0(1D)  O^S)  02  02(xa)  Oat1^)  03 

N N2  NO  N02  N20  NO  3 N205 

H H2  OH  H20  H02  H202 

CO  C02  CH2  CH3  CH4  CHO  HCHO 

CH30  CH302  CH3OOH 

hno2  hno3 


COMPARISON  OF  MEASUREMENTS  AND  CALCULATIONS 

Results  of  some  experimental  measurements  at  higher  latitudes  are  plotted 
on  the  graphs  of  calculated  particle  densities  and  mixing  ratios.  In  many 
cases  referenced  literature  contained  information  on  either  mixing  ratio 
or  number  density,  and  the  other  factor  had  to  be  computed  from  the  avail- 
able data.  When  computation  was  necessary,  the  total  number  density  values 
were  taken  from  tables  in  the  US  Standard  Atmosphere  of  1976. 

For  ozone,  only  the  measurements  by  Randhawa  [5]  on  two  different  days  in 
September  at  Poker  Flat  are  used.  As  can  be  seen  in  Figs,  lc  and  2c,  the 
computational  results  are  very  close  to  the  experimentally  measured  values. 
Although  calculated  values  for  ozone  at  32°  N latitude  [2]  are  not  shown 
here,  the  calculated  values  at  65°  N latitude  were  found  to  be  generally 
higher  than  at  32°  N latitude  in  the  lower  stratosphere  and  lower  in  the 
upper  stratosphere  (25  to  50  km). 

Nitric  oxide  calculated  values  are  shown  in  Figs.  Id  and  2d  with  experi- 
mentally measured  values  by  Ridley  et  al.  [6]  taken  at  59°  N in  July  1974 
and  by  Lowenstein  and  Savage  [7]  in  June  1975.  Comparison  with  Ridley's 
measurements  show  the  calculated  NO  profile  to  be  somewhat  high.  However, 
when  calculated  values  of  NO  are  compared  with  Lowenstein  and  Savage's 
measurements,  agreement  is  quite  good  at  the  altitudes  where  those  mea- 
surements were  made. 

In  Figs,  le  and  2e  calculated  values  of  nitrogen  dioxide  are  fairly  close 
to  measured  values  of  Kerr  and  McElroy  [8],  The  measurements  of  August  17 
were  taken  in  the  evening  and  those  of  August  18  were  recorded  in  the 
morning.  Comparison  with  the  measurements  of  Evans  et  al . [9]  shows  good 
agreement  except  for  a shift  in  the  region  of  maximum  values  and  possibly 
a sharper  decrease  in  N02  concentration  above  30  km. 
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Rate  constants  for  photolytic  reactions  and  for  two-body  reactions  are  given  in  cm3sec_1  and  for  three-body  reactions  units  are  cm6sec 


For  nitric  acid  in  Figs,  lh  and  2h,  the  calculated  values  above  20  km 
agree  quite  well  with  measurements  of  Evans  et  al.  [10];  but  at  20  km  and 
below,  the  measurements  of  Evans  et  al.  and  of  Lazrus  and  Gandrud  [11] 
are  somewhat  lower  than  the  calculated  values.  A shift  in  the  area  of 
peak  concentration  between  measured  and  calculated  values  is  also  observed. 

Calculated  values  of  nitrous  oxide  in  Figs,  li  and  2i  are  shown  to  agree 
with  measurements  taken  by  Schmeltekopf  et  al.  [12]  in  early  spring 
(May  1976)  and  not  quite  as  well  with  measurements  taken  in  late  summer 
(August  1975).  Measurements  in  May  of  1976  were  taken  in  Alaska,  and 
measurements  in  August  of  1975  were  taken  at  Saskatchewan. 

Although  there  is  some  information  in  the  literature  for  experimental 
measurements  of  other  species  in  the  stratosphere  such  as  CH4,  CO,  H20, 

C02,  and  02(1D),  no  measurements  were  plotted  on  the  graphs  of  calculated 
profiles  for  these  species.  However,  differences  in  the  calculated  pro- 
files of  H20,  CO,  and  CH4  between  32°  N [2]  and  65°  N seem  to  agree  with 
comments  on  the  latitudinal  variation  of  these  species  by  Farmer  [13]. 

The  calculated  values  of  water  vapor  are  well  within  the  envelope  presented 
from  measurements  by  Mastenbrook  in  [14]  and  are  extremely  close  to  the 
average  profile  shown  by  Mastenbrook. 

The  authors  do  not  intend  to  slight  measurements  performed  by  experimenters 
other  than  those  mentioned  in  this  report.  Some  measurements  were  not 
included  simply  to  avoid  crowding  the  presented  graphs  to  such  a degree  as 
to  make  them  unreadable.  The  intent  is  merely  to  compare  some  of  the  cal- 
culated results  with  some  observations  to  show  that  the  calculations  are 
somewhat  reasonable  and  to  establish  credibility  in  the  modeling  effort. 
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